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' Abstract We report on results from H band imaging observations of a complete 

^— I , sample of high-luminosity low-redshift QSOs. The luminosity function 

I ■ of QSO hosts is similar in shape to that of normal galaxies, although 

CO ' offset in normalisation by a factor of 10~*. This supports the hypoth- 

, esis that the parent population of quasars is identical to the general 

' population of early-type field galaxies. 

Ph| 1. Introduction 

2 . Identifying the parent population of QSO host galaxies is one of the 

I fundamental problems linking the QSO phenomenon to galaxy evolution 

d ■ in general. While the most widely adopted approach is to compare 

^ i morphological (and increasingly also spectral) properties of normal and 

active galaxies, we investigate here the statistical distribution of galaxy 
luminosities. This short report highlights some results which will be 
I elaborated in more detail in a series of forthcoming papers. We adopt 

h = 0.5 and qo = 0.5 throughout. 

2. H band imaging of a complete QSO sample 

In order to constrain the overall luminosity distribution of QSO hosts, 
the investigated targets must represent a fair sample of the QSO pop- 
ulation. To obtain also the normalisation, the sample has furthermore 
to be complete (in the sense of comprising all QSOs within a given area 
that are brighter than a well-defined flux limit). Our sample has been 
selected from the Hamburg/ESO bright quasar survey (Wisotzki et al. 
2000), with the following additional criteria: Right ascension between 
9^ and 16*^; nuclear absolute magnitudes Mbj < —23; and redshifts 
z < 0.3. The resulting sample consisted of 30 targets distributed over 
2200 deg^, all of which were subsequently observed. Because of the op- 
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Figure 1. Host galaxy versus nuclear absolute magnitudes. Open symbols: disk- 
dominated hosts; filled dots: elliptical hosts. The four radio-loud QSOs are marked 
as circled dots. 



tical selection, the sample is predominantly radio-quiet, containing only 
four bona fide radio-loud quasars. 

Observations were conducted in February 1999 with the ESO NTT 
and its near-infrared camera SOFI. With a seeing between (X'S and C/.'S, 
the host galaxy was clearly detected in all cases, and after PSF sub- 
traction, unique morphological type assignment was possible for most 
objects. A summary of the results: 

■ A disk model is preferred in 11 objects (37%); there is often evi- 
dence for an additional substantial bulge component. 

■ A spheroidal model is preferred in 16 objects (53 %); for the QSOs 
with Mbj < -24, 10 out of 12 (83 %) are spheroidal. 

■ Only three systems (10%) are strongly interacting or irregular. 



■ All hosts are very luminous, Mh ^ M^j with Mh — —24.5 being 
a typical value for the general galaxy population (cf. below). 



The luminosity function of QSO host galaxies 



3 



■ There are no hosts with Mh ^^^ > nuc; thus the McLeod et al. 
(1999) diagonal boundary is confirmed and reinforced (in particu- 
lar, there are no more upper limits). 

■ Except for this boundary, there is no convincing evidence for any 
correlation of Mh ^gai with Mb, nuc (apart from the two extreme out- 
liers 3C 273 and PKS 1302-102, which however are fiat-spectrum 
radio sources and probably beamed). 

3. Bivariate luminosity function 

In order to fully describe a QSO host galaxy luminosity function 
(QHGLF), at least two independent variables are required: Host galaxy 
luminosity Lgai and nuclear luminosity Lnuo leading to a bivariate lumi- 
nosity function $(Lgai, Lnuc)- Cosmological evolution of the QSO popu- 
lation demands <I> to depend additionaly on redshift z. 

Prom the above results, we simplify the expression by tentatively as- 
suming that QSO and host galaxy LFs are formally uncorrelated over 
the luminosity range of the sample under consideration. 



We adopt common parametric forms for both and V'- The former can 
be well approximated, for low redshifts and not too low luminosities, 
by a single power law (Kohler et al. 1997; Wisotzki 2000), the latter is 
usually expressed by a Schechter function. For the evolution we assume 
a simple power-law form of pure number evolution, which is perfectly 
sufficient at low redshifts. The resulting expression for the QHGLF is 



An important additional feature is the McLeod et al. (1999) boundary, 
which basically states that for a given host galaxy mass (and hence lumi- 
nosity), there is a well-defined maximum nuclear luminosity, essentially 
given by the Eddington limit of the central black hole. Although the 
universal validity of this boundary has recently been questioned by Per- 
cival et al. (2000; also these proceedings), it is clearly present in our 
data. We therefore decided to incorporate it by multiplying the above 
expression by a factor ■^(Lgai — a^-^nuc) where TC{ti — i2) = 1 for ti > t2 
and elsewhere is the Heavyside step function, and x is an adjustable 
parameter. For Lgai measured in the H band and Lnuc measured in B, 
we adopt a; ~ 1 (this corresponds to the diagonal dotted line in Fig. 1). 

Numerical values of the QHGLF parameters have been estimated by 
maximum likelihood fitting of the above functional form to the observed 
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Figure 2. Estimates of the QSO host galaxy luminosity function. The points give 
the binned QHGLF with errors and upper limits from Poisson statistics. The lower 
solid curve is the best-fit maximum likelihood estimate. Also shown are some LPs of 
early- type field galaxies for comparison. Upper solid curve: Lin et al. (1999); dashed 
curve: Kochanek et al. (2001). The dotted line is the Lin et al. LP with rescaled 
normalisation to approximately match the QSO data points. 



distribution of sources in {Mb, nuc, Mh^^i^i, z) space. The resulting best- 
fit values are 

1.3 X 10"'^ Mpc~^ for M% Q = -23 
6.5 
-2.8 

-24.6 ± 0.3 
0.5 ±0.5 

We have also determined binned estimates of the QHGLF under the 
assumption that the QSO luminosity function term (^(Lnuo^) is well- 
described by the above parameters. The resulting binned and maximum 
likelihood estimates are shown in Fig. 2. 

The last two parameters are particularly interesting, as they describe 
the shape of the galaxy luminosity function. In Fig. 3 we show confidence 
contours for these parameters and compare them with published values 
determined for the general field galaxy population (adapted to the H 
band by assuming typical galaxy colours). The characteristic luminosity 
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Figure 3. Likelihood contours for the HGLF parameters: The solid line corresponds 
to 68 %, the dashed line to 90 % confidence. Comparison values for the field galaxy LF 

arc given by symbol markers. Open circle: Lin ct al. (1999); filled circle: Gardner ct 
al. (1997); open square: Kochanek et al. (2001); filled square: Loveday et al. (2000). 



gal completely within the range found for field galaxies, while the 
faint-end slope (5 is only slightly flatter. The latter makes the QHGLF 
formally inconsistent with that of field galaxies, but Fig. 2 shows that 
the discrepancy even at the faint end is not large: The dotted line rep- 
resents the Lin et al. (1999) LF of early-type field galaxies, rescaled to 
approximately match the normalisation of the QHGLF data points. The 
overall agreement is surprisingly good, and the slight mismatch at the 
faint end might even be explained by our limitation to only luminous 
QSOs. 

To summarise, our results present no evidence that QSO host galaxies 
show a luminosity distribution widely different from that of inactive field 
galaxies. This lends considerable support to the hypothesis that the 
parent population of quasars consists of ordinary early-type galaxies, 
with very little if any bias towards higher luminosities or higher degree 
of morphological irregularities. 
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4. Implications for the QSO duty cycle 

Dividing the space density of quasar hosts by that of the parent galax- 
ies gives the ratio of active to inactive galaxies as a function of galaxy 
luminosity. This ratio has been interpreted in the past as the 'probabil- 
ity' of finding a QSO in a galaxy. By assuming that all major galaxies 
have massive black holes in their centres and using the argument of 'en- 
semble average equals time average', this is furthermore the time fraction 
that a typical galaxy spends in a 'quasar state' (also called the QSO duty 
cycle 6). 

Prom our determination of the QHGLF, we found that the parent 
population can be probably identified with normal early-type field galax- 
ies. Under these premises, we can estimate S directly from the rescaling 
factor used in Fig. 2 to match the field galaxy LF to the QSO data 
points, giving S ~ 10"^. Because of the similarity between the lumi- 
nosity functions, there is virtually no trend of S with luminosity. We 
can confidently exclude that the number density ratio between quasar 
hosts and luminous field galaxies approaches unity (implying 5 ~ 1) as 
suggested recently by Hamilton et al. (2001). 

Note that the above value for S refers to zero redshift, since because 
of the explicit evolution term for the QSO space density in the analytic 
expression for the QHGLF normalisation is taken at 2; = 0. In this 
simple picture, QSO evolution would be interpreted as an increase in 
the duty cycle towards higher redshift with a rate of 5 ex (1 -|- z)^'^. 

With < 2; < 0.3 for our sample, this corresponds to a time span of 
T ~ 5Gyrs. Multiplying this with the redshift-avcragcd value of S{z), 
we can convert the duty cycle into a time scale of < At >~ 10^ yrs, which 
is a robust lower limit to the total time a typical low-redshift quasar is 
switched on, and probably not too different from the actual time scale 
for quasar-type activity. 
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